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Abstract 

We introduce a general method which allows reconstruction of electronic band structure of 
nanocrystals from ordinary real-space electronic structure calculations. A comprehensive study 
of band structure of a realistic nanocrystal is given including full geometric and electronic relax- 
ation with the surface passivating groups. In particular, we combine this method with large scale 
density functional theory calculations to obtain insight into the luminescence properties of silicon 
nanocrystals of up to 3 nm in size depending on the surface passivation and geometric distortion. 
We conclude that the band structure concept is applicable to silicon nanocrystals with diameter 
larger than !^ 2 nm with certain limitations. We also show how perturbations due to polarized 
surface groups or geometric distortion can lead to considerable moderation of momentum space 
selection rules. 
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I. INTRODUCTION 



Crystalline nanostructures are often viewed as "artificial atoms" or "zero-dimensional 
systems" despite being the size of a smaller protein. This is because when single nanocrystals 
(NCs) are experimentally probed, very sharp optical radiative transitions, implying the 
existence of discrete energy levels rather than energy bands, are observed. [1] On the other 
hand, nanostructures are known to retain some of the band-structure-related properties of 
their bulk counterparts, e.g. direct-bandgap semiconductor NCs such as CdSe are excellent 
light emitters. 

This brings about an apparent dichotomy and raises the question if, and under what 
circumstances, the band structure concept can be applied to NCs and, moreover, if the 
band structure description can overlap with the description based on discrete energy levels. 
In this article, we show that band structure can be rigorously defined in nanostructures, 
taking into account several alterations with respect to bulk materials. 

The concept of band structure in NCs is very important, because it allows one to transfer 
the already thoroughly-studied electronic properties and the whole framework of solid-state 
physics of bulk materials to the nanoscale. It might not play such a crucial role in direct- 
bandgap materials, which have already been successfully theoretically described using models 
completely disregarding the band structure, [2] but it starts to be of particular importance 
for the case of indirect-bandgap semiconductors. 

In indirect-bandgap semiconductors, the comprehension of the electronic band structure 
behavior allows for efficient band engineering. A prime example illustrating the usefulness of 
the band-structure approach is germanium, an originally indirect-bandgap material success- 
fully transformed to a direct-bandgap one by heavy doping and tensile strain. [3] This exper- 
imental realization of direct-gap bulk germanium was based on theoretical band-structure 
calculations. [4] What is important, however, is that the same concept was experimentally re- 
alized also in germanium NCs. [5] This example confirms that the parallels between bulk and 
nanocrystalline materials do exist and that they can be beneficial in material engineering as 
well as for the intuitive understanding of the material's behavior. 

It is evident from the discussion above that in direct-bandgap materials only the states 
close to the valence and conduction band edges (situated at the F point) will be those govern- 
ing the luminescence behavior of the material. However, in indirect-bandgap semiconductors 
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band structure can be influenced in a more complex way, making it more diflicult to decipher 
the resulting more complex behavior. This is why we chose to analyze the band structure 
of nanocrystalline silicon, as silicon is the most wide-spread indirect semiconductor. 

Silicon as material found wide applications in electronics and it is still a subject of intense 
research. Nevertheless, the indirect band gap nature of the band structure of bulk silicon 
has always been the major obstacle for its employment in light-emitting devices since mo- 
mentum conservation requires additional momentum transfer mechanisms involved in the 
light emission processes. The situation changed dramatically in the last two decades due to 
the emergence of the possibility of preparing Si-based structures with nanometer size, where 
quantum effects begin to play a dominant role. In particular, large effort has been devoted to 
the study of the optical properties of Si NCs in the last years, with a perspective of potential 
for real-life applications such as e.g. light emitting diodes, next-generation solar cells and 
biomedical devices. [6-9] The discovery of efficient visible photoluminescence[10] (PL) and 
optical gain[ll] from sihcon NCs has demonstrated the possibility of partially overcoming 
the limitations of the indirect band gap of silicon by exploiting the quantum phenomena at 
the nanoscale. Despite the large amount of papers published on this subject, there are still 
many aspects which are not fully understood and are subject of intense dispute. 

As already mentioned, the theoretical concept of optical properties of Si NCs is often 
discussed in the framework of the band structure picture of the bulk material. However, 
the flnite size of a system measuring only a few nanometers makes the justiflcation of this 
approach questionable. Therefore, the validity of the band structure concept and, if need 
be, its character in Si NCs of a given size and surface passivation have to be analyzed. 

Traditional electronic structure calculations of Si NCs, as a flnite system, provide only 
real space molecular orbitals (MOs), where crystal momentum and thus band structure E(k) 
is not directly accessible. What's more, in NCs with typical sizes of a few nanometers, the 
surface-to-volume ratio substantially increases and thus the interface between the NCs and 
their environment plays a crucial role in tailoring their optical properties. Consequently, the 
optoelectronic properties of Si NCs are very sensitive to their surface passivation, symmetry 
and applied strain, which makes the application of simple (bulk like) models questionable. 

To the best of our knowledge, only a few attempts[12-15] have been made to perform 
the projection of states of flnite systems from real to reciprocal space and vice versa so far. 
Therefore, a robust method which allows the band structure mapping from fully relaxed 
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density functional theory (DFT) calculations of realistic Si NCs is required. 

From the perspective of theoretical simulations, many approaches, ranging from ab 
initio[16-20] to parametrized semiempirical methods [12, 21-23] , have been adopted to 
investigate the optical and electronic properties of Si NCs. Ab initio methods provide a very 
accurate description of electronic states of fully relaxed Si NCs, but simulations of realistic 
systems with nanometer size are impossible due to excessive computational demand. On 
the other hand, semiempirical methods [12] might address Si NCs of realistic size consisting 
of tens of thousands of atoms, but with limited transferability. From this point of view, the 
application of fast local orbital DFT codes [24] devised with the aim of computational effi- 
ciency but still providing desired precision in description of ground state electronic structure 
seems to be the optimal choice. 

In this paper, we will introduce a general method which allows the reconstruction of band 
structure from MOs obtained e.g. from ab initio calculations. For this reason, we dare to 
use the terms highest occupied molecular orbitals (HOMO) and lowest unoccupied molec- 
ular orbital (LUMO) and valence-band maximum/conduction-band minimum as synonyms 
throughout this article. We will demonstrate the method's capability to fold up the band 
structure of finite systems on a simple ID atomic chain consisting of a few H atoms. Using 
this example, we will discuss the main characteristic features of the band structure of finite 
size systems. Next, we will apply the procedure to analyze the electronic structure of free- 
standing Si NCs obtained from DFT calculations. To achieve this objective, we will adopt 
fast local orbital DFT code Fireball [24? ], which allows us to perform fully relaxed total 
energy calculation of different Si NCs consisting of up to thousand of atoms in a feasible 
way. We will show how the electronic band structure of Si NCs is affected by their surface 
passivation, symmetry and size. 

II. METHODS 

A. DFT calculations 

All computations were carried out using the local orbital density functional theory (DFT) 
code FIREBALL [24? , 25] within the local-density approximation (LDA) for the exchange- 
correlation functional. Valence electrons have been described by optimized [25] numerical 



4 



atomic-like orbitals having the following cutoff radii (in a.u.): Rc(s, s*) = 4.0 for H; Rc(s) = 
4.5, Rc(p) = 4.5, and Rc(d) = 5.4 for C; Rc(s) = 4.8, Rc(p) = 5.4, and Rc(d) = 5.2 for Si, 
and Rc(s) = 3.5, Rc(p) = 4.0, and Rc(d) = 5.0 for O, respectively. The correctness of the 
basis set was checked to reproduce the band structure of bulk silicon within LDA accuracy. 

Nanocrystals are represented by cluster models with 3 different core sizes consisting of 
68, 232 and 538 Si atoms, respectively. These 3 models represent Si NCs with the diameters 
of 1.5, 2.0 and 2.5 nm Fig.l. 

Si core of each of the three models is terminated either with polar (-0H, simulating an 
oxidized layer) or non-polar (-H, -CH3) passivating groups. Atomic structure of the Si^ss 
core was cut out from relaxed bulk silicon lattice in such a way so as to minimize the number 
of unsaturated bonds and to reflect the lattice symmetry with well-defined (111) and (100) 
faces ( Fig.l ). Subsequently, the smaller Si232 and Sies cores were derived from Si^^s by 
removing the topmost atomic layer, in order to get a smaller analog of the same symmetry 
and surface faces. Then, passivating groups of the three types were attached to the Si core, 
leaving all surface Si atoms fully saturated. 

All Si NCs models were fully optimized, allowing the relaxation of all atoms. Total 
energy calculations were performed as a cluster calculation ( k=0 ) and the convergence was 
achieved when the residual total energy of 0.0001 eV and the maximal force of 0.05 eV/ A 
were reached. 

B. Momentum space projection of molecular orbitals 

The band structure theory of solids is well established and it has been successfully applied 
to real materials to explain their physical and material properties, e.g. electrical resistivity 
and optical absorption. A corner-stone of the band structure theory is the so-called Bloch's 
theorem, [26] which had initiated the epoch of modern solid state physics. Electronic states 
in an infinite periodic system are described via the Schroedinger equation: 

[V' + y(r)]^^(r) =e*(r), (1) 

where e is the energy of eigenstate *(r) and potential V{r) is periodic V{r) = y(r + R), 
R being the translational lattice vector. Bloch's theorem postulates that the solution of 
the Schroedinger equation, ^^k(r), can be written as a product of a real mother- function 
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FIG. 1. Geometry of model nanocrystals. First two panels show 9 different models of NCs with 3 
core sizes (68, 232 and 538 Si atoms of diameter 1.5, 2.0 and 2.5 nm) with either polar oxidized 
(-0H) or non-polar (-H, -CH3) passivation. The last panel shows the geometry of Sisas-H from 
important crystallographic directions (100), (111) and perpendicularly to (111). 

'^n,k(r), which is also periodic t^ri,k(r) = '?in,k(r + R), and a Bloch plane- wave e^*^-^: 

*n,k(r)=^ln,k(r).e^'^■^ (2) 

The eigenstate ^^ri,k(r) varies continuously with the wave vector k and forms an energy band 
e^^k identified by the band index n. 

How does the band-structure picture change when the size of the system is reduced? Or, 
in other words, is the concept of energy bands also valid for nanoscopic systems? As we will 
show in the following, the finite size of a system has two important consequences: (i) wave 
vector k becomes discrete and (ii) wave functions are delocalized in momentum space [27]. 

The first statement (i) directly follows from the construction of reciprocal space. For 
example, let us assume a 1-D mono-atomic chain of atoms with the lattice constant a. 

6 



There are N different wave vectors k separated by 

" a.iV " L ' ^"^^ 

where L — N.a is the length of the chain. Obviously, if N becomes small, the separation 
between the wave vectors Ak^^^ becomes larger. 

Secondly, according to the Heisenberg uncertainty principle the crystal momentum be- 
comes delocalized as follows: 

AxAp > h. (4) 

Substituting the crystal momentum P — and the size of the nanocrystal Ax = L into 
the equation above, we obtain the relation for the delocalization of wave vector due to the 
finite size of the system 

Ak('^ > ^. (5) 

The fact that the separation of discretized wave vectors Afc^^^ and their delocahzation Afc^^^ 
are of the same order calls for a more rigorous discussion of the band structure of finite 
systems, which is provided in the following text. We will use natural units {h = 1), using k 
as a synonym for momentum. To simulate the finite size of systems (e.g. a nanocrystal) we 
introduce the so-called window function w{r) which restricts the wave function ^n,k(r) to 
the space occupied by the system. Then, Eq.2 is modified accordingly: 

*„,k(r)=w;(r).M„,k(r).e^'^-'-. (6) 

Fig.2 depicts the characteristics of the individual terms in Eq.6 in a 1-D case. In the simplest 
approximation, one can define the window function w{r) as a stepwise function 

. 1 : |r| < L 

w{x) = { (7) 
: |r| > L, 

which vanishes outside the nanocrystal (see Fig. 2a). Then, the resulting real-space wave 
function ^n,k(r) consists of the modulation of the mother function by the Bloch plane wave 
and the window function w{r) (see Fig.2d). In real systems, the window function can have 
a more complicated shape dictated by the electronic structure of a particular system (see 
green dashed line in Fig.2a). 
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FIG. 2. Schematic illustration of Eq.6. (a) Window function constraining the wave function to 
a finite area of space. In general, it is the envelope of a MO and it can have a complicated 
shape (green dashed line). However, for simplicity, we expect the rectangular shape of Eq.7 (gray 
rectangle), (b) Mother function of an infinite system composed of a linear combination of atomic 
orbitals. (c) Infinite Bloch plane wave, (d) Final wave function of a nanocrystal as a product of 
(a),(b) and (c) constrained in a finite area of space. 

Let us consider a problem inverse to the computation of band structure from Bloch's the- 
orem (Fig. 3a). We know the real-space wave function (/>i(r) of i-th eigenstate (e.g. a Kohn- 
Sham MO) obtained from an aperiodic DFT calculation of a nanoscopic system) and we 
would like to assign the corresponding crystal momentum k. Each MO of energy can 
cross one or more bands e^^k (Fig. 3b). Then, we can assume the MO (/>i(r) to be a linear 
combination of A bands ^^^^k(r): 



,(r) = 5^4^,,,k«(r), 



(8) 



a=0 



where, for a particular 0i(r), neither the indexes n^^k^ nor the coefficients da are explicitly 




FIG. 3. Schematic illustration of Eq.8. (a) Usual approach to obtaining the band structure of an 
infinite crystal using Bloch's theorem. For a given k-vector, the energy spectrum of e^^k is found, 
(b) Inverse approach in a finite system. For a given energy state of orbital 0^, the corresponding 
k- vectors need to be found. The orbital is a linear combination of all Bloch-like wave functions 
^n,k(r) corresponding to the n, k in which the orbital energy crosses the energy bands en(k). 

known. 

One option how to extract k- vectors of (/>i(r) lies in the determination of the corresponding 
^na,ka(^) t^na,ka(^) E consecutive extraction of the Bloch factor e^*^-^ from Eq.6. This 
is possible, for example, by projecting (/)i{r) on ^^^j^^^(r) obtained from bulk calculation 
[28]. Nevertheless, wave functions ^^ri,k(r) in a finite system may differ considerably from the 
bulk one ^^^^^^(r) (e.g. if w{r) is more complicated) and a large number of states have to 
be considered for a large nanocrystal. These complications hamper a robust implementation 
of this method to a computational code. 

Here we opt for another method, which transforms MO 0i(r) from real into reciprocal 
space using Fourier transform. A similar approach was successfully applied to analyze an- 
gular resolved photoemission (ARPES) spectra of oligomeric organic molecules [13]. In our 
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scheme, selected ^^(r) expanded on a real space grid is projected to momentum space either 
using discrete 3-D Fast Fourier Transform (FFT) or by the projection onto a set of plane 
waves e . 

0,(k') = (0,(r)|e^'^'"^), (9) 

where k' is an arbitrary wave vector. 

For practical use, we plot the so-called momentum density y5^(k'), which is real, instead of 
the complex Fourier transform (/)i(k'). The momentum density y5^(k') can be written, using 
Eq.9, as follows: 

p,(k') = (10) 

It also contains the information about the delocalization of the MO in momentum space and 
it can be plotted as a function of wave vector along a selected line in momentum space (see 
Fig.4d). 

The first approach (FFT) provides a picture of the 3-D structure and symmetries of 
the particular state in reciprocal space (see Fig.4e). Nevertheless, its resolution is limited 
by the size of the real-space grid on which MOs are expanded. The second approach (the 
projection on set of plane waves) is more suitable to stick with the traditional 1-D band 
structure representation plotted along the lines connecting the high-symmetry points in k- 
space. In this case, we let k' sample the given high-symmetry line in k-space with much 
higher resolution. Although this approach is straightforward, a rigorous analysis of the 
structure of the resulting momentum space distribution is required. For simplicity, let us 
assume a one-to-one correspondence between a particular MO and a band wave function 
(/)i{r) = ^/n^\^{r). Note that the generalization of the following discussion 

if 0i(r) is a hnear combination of several *na,ka(i*) (Eq.8) is straightforward due to the 
linearity of Fourier transform. Fourier transform of a wave function ^n,k(r) given by Eq.6 
can be expressed as a convolution of the three terms: 

*,,k(kO = wik')^uik')^5ik' - k). (11) 

It is well worth analyzing in detail the process of convolution and the character of each 
term in Eq.ll. To make our discussion more illustrative, we restrict ourselves to 1-D case. 
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FIG. 4. Schematic illustration of convolution of momentum space representations of the three wave 
function components, (a) Rectangular window function transformed to sinc(Lk'). (b) Bulk-like 
mother function transformed to a series of 5-functions in the centers of Brillouin zones, (c) Bloch 
wave transformed to a 5- function inside the first Brillouin zone shifted out of the center, (d) 
Momentum space representation of the wave function created by the convolution of (a),(b) and (c) 
in a 1-D case, (e) Illustration of a more-D case: a 2D view of momentum space density of MOs 
in a silicon nanocrystal. Note the side artifacts spaced by multiples of K representing the higher 
Fourier components in the expansion of mother function. 

Fig. 4 represents schematically the process of convolution and the character of each term 
separately. 

We start from the trivial third term. 5{k' — k) is a Fourier transform of Bloch wave e^*^-^, 
where k is from the first Brillouin zone (see Fig. 4c) and attains some of the discrete values 
separated by A/c^^^ = 27r/L as postulated in Eq.3. 

The first term in w{\<:') in Eq.ll is a Fourier image of the window function w{y). This 
term causes the delocalization of wave function in k-space. Considering the window function 
w{y) as the step function of length L (see Eq.7), its Fourier transform equals w{k') = 
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sinc(Lk'). Using Eq.5, it can be expressed also as sinc(27rkYAfc(^)), which brings us back to 
the Heisenberg principle. 

The square |^^;(k')p introduces to the momentum space density Pi(k^) a blur with a 
Lorentzian envelope (see Eq.4a) of slow asymptotic decay l/(Lk')^. This has an important 
implication for the momentum selection rules for the transitions between states in finite 
systems: two wave functions centered around different k can have non-negligible overlap in 
k-space even if their separation in k is fairly large. 

We would like to stress that the decay in the momentum space projection depends strongly 
on the actual shape of w{r). Therefore, this approach is superior to a simple estimation of 
momentum uncertainty based on the Heisenberg principle, which can provide only the width 
of the peak around k, but says nothing about its decay. 

Finally, the second term in Eq.ll u(k.^) corresponds to the Fourier transform of the 
mother function t^n,k(r), which can be expanded in a discrete Fourier series of waves using 
its periodicity: 

oo 

Mn,k(r) = 5]c^e^("^•^)•^ (12) 

where K denotes the reciprocal lattice vector and m is an integer index addressing different 
reciprocal unit cells. This expression is transformed to reciprocal space as follows: 

oo oo 

<k(k') = {J2 c^e^^™-^)-"^|e^'^'-'-) = J2 - m.K), (13) 

where S(k^ — m.K) is situated in the center of the m-th reciprocal unit cell (i.e. not in the 
first Brillouin zone). 

Typically, t^n,k(r) is relatively smooth in the unit cell and so the Fourier expansion coeffi- 
cients Cm decay with increasing frequency (Fig.4b). In reality, just one expansion coefficient 
is dominant. The dominant Fourier component is determined by the nodal structure of 
'^n,k(r), where nodes are introduced either by anti-bonding character of the mother function 
or by the contribution of the higher angular momentum atomic orbitals (such as a p-orbital). 
If '^n,k(r) has no node inside a unit cell (for example bonding state of two s-orbitals), the 
first coefficient cq situated in the center of the first Brillouin zone (k' = 0) is the dominant 
term. Similarly, if Un,\^{r) has m nodes, the m-th Fourier expansion coefl&cient Cm dominates. 
This means that most of the momentum space density is situated in the m-th reciprocal unit 
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FIG. 5. 1-D cross-section of momentum-space projected MOs of a chain of 8 hydrogen molecules, 
(a) Atomic structure of molecular hydrogen chain with bond lengths of 0.7 Aelectronically coupled 
through the intermolecular distance of 1.3 A. (b) Real space projection of MOs 0i(r) of the system. 
First 8 orbitals are bonding and occupied, the remaining 8 are anti-bonding and unoccupied. 
Positive sign of wave function is colored in blue, the negative sign in red. (c) Band structure of 
the finite chain as assembled from momentum densities pi(k^) plotted at particular orbital energy 
ei by color scale. Parabolic band dispersion can be clearly seen following the maxima of densities 
(yellow-red-black). For some of the orbitals z = 1, 8, 9, 13, the actual shape of function Pi(k') is also 
plotted to illustrate the correspondence to Fig. 4. The main peaks corresponding to the Fourier 
components in the expansion of t^ri,k(kO are denoted by arrows. 

cell. Consequently, the convolution of Un,\i{^^) with the Bloch term 6(k^ — k) forms a new 
band located in the m-th reciprocal unit cell (see the anti-bonding band in Fig. 5c). Thus, 
our approach provides the reconstruction of the so-called unfolded band structure, where 
each band is situated in a different reciprocal unit cell. Note also that in a 3-D case the 
index m is a vector, and nodes in each dimension should be considered independently. 
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Now, we will illustrate the projection method on a simple case. Let us analyze the 
electronic structure of a finite 1-D chain consisting of 8 coupled hydrogen molecules (see 
Fig. 5), where the unit cell consists of a single hydrogen molecule. Upon projecting real 
space MOs (Fig. 5b) into k-space two distinct bands isolated by a band gap appear (see 
Fig. 5c). The lower band (m = 0) composed of bonding orbitals with no node inside the unit 
cell has the dominant Fourier components located in the first Brillouin zone. The higher 
anti-bonding band {m — 1) has the dominant Fourier components located in the second 
reciprocal unit cell. 

Apart from these dominant Fourier components, which characterize the band structure, 
there are also smaller satellite peaks representing the other components in the Fourier 
expansion of Un,k(k.')- In the special case, when Bloch k = 0, the peaks are located directly 
in centers of m-th reciprocal unit cell and correspond directly to components Cm in Fourier 
expansion of 'S^^k(k') (Eq.l3). This is illustrated for MO z = 16 by peaks denoted Co,Ci,C2 
in Fig. 5c. In general case (k 7^ ) these peaks are split and shifted by ±k as denoted by 
K + k and K — k in Fig.5c. 

To illustrate the delocalization of the momentum vector k due to the finite size of the 
system, represented by w{r)^ we analyzed the lowest band ( i = 1 ) of several different 
1-D hydrogen chains with lengths from 2 to 16 molecules. In Fig.6a the peak width of the 
depicted momentum density decreases proportionally to the number of molecules (unit cells) 
in the chain according to Eq.5. In this particular case, the window function w{r) is very 
close to a rectangular step function. Hence, the shape of the momentum space density p(k^) 
is very similar to |sinc(Lk^)p. The frequency of this sinc-like function can be also deduced 
from the number of nodes per reciprocal unit cell, which is proportional to the increasing 
chain length L. 

Let us summarize the main conclusions of this chapter. We discussed two main implica- 
tions of the finite size of a nanocrystal on its band structure: (i) the discretization and (ii) 
the delocalization of wave vector k. We introduced a robust method of projection of MOs 
to k-space, which provides band structure of a finite system, and discussed in detail the 
analysis of the resulting momentum densities p^(k'). We demonstrated that the important 
characteristics (i.e. Bloch k-vector, its delocalization and the discrimination of independent 
bands) can be extracted from the resulting momentum densities, despite the convolution in 
Eq.ll. 
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FIG. 6. Dependence of the delocalization of momentum density Pi(k') of MOs in a chain of hydrogen 
molecules on the chain length, (a) Momentum space delocalization of the lowest MO {i — 1) in 
a chain composed of 2 (black), 4 (red), 8 (green) and 16 (blue) hydrogen molecules, (b) Band 
structure of a chain composed of 2 hydrogen molecules, (c) Band structure of a chain composed 
of 16 hydrogen molecules. 

C. Band structure of bulk silicon and nanocrystal 

Before we advance from the 1-D hydrogen chain to the band structure of a silicon NC, 
it is worth discussing in detail the reciprocal space structure of bulk silicon. While in solid 
state physics it is common to describe band structure in the compact picture (or reduced 
zone scheme), where all bands are folded into the first Brillouin zone, for our purposes it 
is more suitable to use the so-called unfolded picture. In this approach, bands located in 
higher reciprocal unit cells are considered independently. 

Silicon crystallizes in a diamond lattice consisting of two interpenetrating face centered 
cubic (FCC) Bravais lattices. In real space, the Wigner-Seitz cells of FCC build up a rhombic 
icosahedral honeycomb, which converts to a truncated octahedral honeycomb in reciprocal 
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FIG. 7. 3D structure of silicon reciprocal lattice, (a) Structure of the first Brillouin zone of silicon 
in the shape of a truncated octahedron with high symmetry directions shown, (b) Positions of 
points Ta (blue), (red) and X55 (green) around the first Brillouin zone, (c) The points F^, F5 
and X55 in the context of the sub-lattices a(blue) and b (red), (d) Reciprocal unit cells adjacent 
to the first Brillouin zone, distinguished to the sub-lattice a (blue) and b (red), (e) 3D view of the 
localization of maximal k-space projected density of HOMOs (green) and LUMOs (purple) for a 
Sies-H NC. Points F^, F5 and X55 marked by spheres to highlight the correspondence. 

space. As shown in Fig.7d, the reciprocal lattice is composed of two cubic sub-lattices a 
(blue) and b (red) with a truncated octahedral shape. The sub-lattice a is centered in Fq 
of the first Brillouin zone and the sub-lattice b is shifted by a vector (1,1,1) with respect to 
the sub-lattice a. Thus, the sub-lattice b is situated in the cube vert exes and the sub-lattice 
a in the cube centers of body centered cubic lattice (BCC), as shown in Fig. 7b. 

The unit cell located at Fq is the first Brillouin zone. There are in total 14 reciprocal 
unit cells adjacent to the first Brillouin zone. Six of them belong to the sub-lattice a, being 
centered at the points F^ with coordinates V (±2,0,0). Another 8 unit cells belongs to the 
sub-lattice b centered at F^ with coordinates V (±1,±1,±1). Here, V means all permutations 
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of axes and signs (e.g. (±2,0,0), (0,±2,0) and (0,0,±2) ). 

All these 14 reciprocal unit cells are of crucial importance for us, because both the valence 
and conduction bands are situated there (in unfolded picture). In particular, the valence 
band maximum is located in the points and F^. The absolute conduction band minimum 
in bulk silicon is located near the points of coordinates V (±1,±1,0), which are halfway 
between two F5. The positions of the points F^, F5 and are schematically depicted 
in Fig.Tb. The localization of maximal k-space projected density obtained from HOMOs 
(purble) and LUMOs (green) for a Sieg-H NCis shown in Fig.Te. From this figure, it is 
evident that HOMOs are localized at F^, F5 and LUMOs near the points. To make 
the context more clear, in Fig. 7c we also depict the positions of the important points F^, 
F5 and in the context of the sub-lattices a and b. Because of convolution with Fourier 
expansion of mother function ( Eq.l3 or more illustrative Fig.4b,d ) both the valence and 
the conduction bands are situated outside of the first Brillouin zone in unfolded picture. 
Therefore, the first Brillouin zone does not play any significant role in optical transitions 
and it is depicted in gray in Fig. 7b, c. 

In the case of bulk silicon, all these reciprocal unit cells around F^ and F5 are equivalent 
due to lattice symmetry. Therefore also the frontier orbitals (HOMOs and LUMOs), which 
are composed of Bloch states located in F^, F5 and are energetically degenerate. More 
precisely, HOMOs (resp. LUMOs) are arbitrary hnear combinations of Bloch states from 
different F^, F5 (resp. X55) as was discussed in previous chapter ( Eq.8 and Fig.3 ). 

However, this is no longer true in NCs, where the symmetry is broken due to different 
lengths in each direction or due to anisotropic perturbations (e.g. mechanical strain) induced 
by surface passivation. This is another reason why it is necessary to describe the band 
structure of a nanocrystal in the unfolded picture with all reciprocal unit cells independent. 

In other words, the anisotropy of NCs removes the degeneracy of MOs. We can observe 
this in our model NCs even though the shape of cores (Sies, Si232, S^ss) respects the cubic 
symmetry of silicon lattice. In the case of -H passivation, the degeneracy is almost preserved, 
because no anisotropic strain or electrostatic field is induced. The LUMO state is 12-fold 
degenerate, representing the 12 points 

However, due to the shape of a Sies-H nanocrystal — slightly shorter in the (111) direction 
(Fig.Sa), MOs are split into two classes of slightly different energies. Electronic states with 
k- vector oriented closer to (111) have energy of 2.17 eV, while states with k- vector oriented 
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FIG. 8. Degeneracy loss due to the broken symmetry of a Sies-H nanocrystal. (a) Top view of a 
Sies-H nanocrystal along the Ill-direction, (b) Side view of a Sies-H nanocrystal perpendicularly 
to the Ill-direction, (c) Two classes of energy split MOs of Sies-H due to anisotropic confinement. 
There are 6 degenerate states in the directions closer to Ill-direction (polar, blue) with the energy 
of 2.17 eV, and other 6 degenerate states more perpendicular to the Ill-direction with the energy 
of —2.03 eV (equatorial, red). Note that states k and — k are always degenerate due to time reversal 
symmetry and the realness of wave function in a finite system. 

more perpendicularly to the (111) direction have energy of —2.03 eV (see Fig. 8c). 

Finally, we make a comment on the importance of MO symmetry for optical transitions. 
Since MOs are real functions, their Fourier transforms are symmetric for k and — k. Con- 
sequently, we can identify two distinct classes of MOs according to the character of their 
symmetry: (a) if a MO is an even function in real space, its k-space representation becomes 
real (cosine part); (b) if a MO is an odd function, the k-space representation is purely imag- 
inary (sine part). According to selection rules, optical transitions are strongest between two 
states localized around the same point in k-space with complementary symmetry character, 
i.e. sine to cosine or vice versa. 
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FIG. 9. Projected density of states around gap split into the contributions of surface (red) and 
core atoms (blue), for the non-polar passivation (a) Sisas-H and (b) Si538-CH3; frontier orbitals 
are localized mostly inside the core. However, for an oxidized NC (c) Sisas-OH, these are mostly 
surface states, (d) shows variation of energy HOMO-LUMO gap with core diameter for all three 
types of passivation 



III. RESULTS 



In this section, we will employ the method introduced in the previous section to analyze 
the band structure of Si NCs as a function of their diameter and surface passivation of Si 
NCs of nanometer size. Namely, we will analyze an impact of surface passivation including 
both polar (^OH) and non-polar (-H,-CH3) groups on their atomic and electronic structure. 



A. SiNC:H 



First, we consider Si NCs passivated by hydrogen groups (^H). Our calculations of fully 
optimized structure show that the atomic relaxation of Si atoms in the core region is negli- 
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FIG. 10. Spectrum of electronic states (the black, green and red bars) with their k-space projection 
from r to X (color scale legend is on the right) for (a) Sisas-H, (b) Sisag-CHa and (c) Sisas-OH. 
The real space projections of selected frontier orbitals are plotted as well in order to show their 
overall shape, symmetry and localization. Note: colorscale was limited to the same maximal value 
in order to enhance the contrast and to compare the absolute k-space densities between different 
passivations. 

gible. Therefore, the Si NC core has very similar atomic structure as bulk Si. We attribute 
this effect to a negligible charge transfer between the Si core and the hydrogen terminat- 
ing groups due to fairly close electronegativities of both elements (1.9 for Si, 2.2 for H, 
respectively) and negligible mechanical stress resulting from a small volume of the hydrogen 
terminating groups. 

Calculated projected density of states (PDOS) for a Si^ss NC is shown in Fig. 9a. Surface 
states corresponding to Si-H bonds are localized far from the band gap edges. Therefore, 
highest occupied and lowest unoccupied electronic states are localized mainly in the core 
region. The electronic states in the Si NC core are strongly affected by the quantum confine- 
ment effect. Consequently, the band gap of a Si NC decreases with increasing NC diameter 
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FIG. 11. Comparisons of 1-D cross-sections of projected k-space density in -H passivated NCs of 
different sizes for (a) Sies-H, (b) Si232-H and (c) Sisas-H. Gradual emergence of energy bands can 
be seen with increasing NCsize. (d) Band structure of bulk Si for comparison. Note: color scale 
was normalized by maximal density independently for each NC size. 

(see Fig.9d), in good agreement with previous calculations [16, 18-20]. 

Fig. 10a displays real space wave functions corresponding to individual electronic states 
on band gap edges. We see that, in case of hydrogen passivation, both HOMOs and LUMOs 
are spread almost homogeneously over the whole Si^ss core. Applying our transformation 
method, we can convert electronic states from real to reciprocal space to obtain the electronic 
band-like picture of a Si NC with a given diameter. Fig. 11 represents projected k-space 
density along the T — X direction. We can see the gradual emergence of the band-like 
structure with the increase of NC size. While for the smallest Sies NC the band structure 
can be hardly recognized, in the case of the Si^ss NC the dispersion of electronic states near 
the band gap mimics well the bulk band structure of both the conduction and the valence 
bands (compare Fig.llc,d). Therefore, we can estimate ^2 nm to be a phenomenological 
limit where it makes sense to speak about electronic band structure and indirect band-gap 
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FIG. 12. Delocalization of k-space projection of LUMO in -H passivated NCs depending on the 
nanocrystal diameter. The sinc-hke shape can be clearly seen, which means an approximately 
rectangular envelope of the MO in real space. 

in silicon NCs. 

There are two important differences between bulk band structure and projected band 
structure of a Sisas-H NC that has to be noticed: (i) presence of discontinuities (mini-gaps) 
in energy dispersion and (ii) delocalization (blurring) of electronic states in k-space. 

The presence of the mini-gaps within the bands is induced by the finite size of the Si 
NC and the preserved symmetry of the atomic Si-core structure. The latter evokes the 
degeneracy of the electronic states with k-vector in the equivalent crystal lattice directions 
in a similar manner as shown in Fig. 8 for Sieg-H. We should note that the presence of 
mini-gaps larger the 63 meV can have important implications for the relaxation process of 
hot excited electrons from F to X. Usually, the maximal vibrational energy of phonons in 
bulk Si does not exceed 63 meV [29] . Therefore hot electrons might be unable to reach the 
conduction band minimum near the X-points once they meet the mini-gap along their path 
from F to X. 
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FIG. 13. Map of Hatree potential inside Si^ss NCs with different passivations. For non-polar 
passivating groups (a) -H and (b) -CH3, the potential is regular and corresponds to the charge 
distribution in the crystal lattice. However, for polar -OH passivation (c), considerable local 
electrostatic fields can be seen near the surface. Note: Potential of neutral atoms was substracted 
from the full self-consistent potential to make the features easily visible. 

The second difference between the Si-bulk band structure and a hydrogen passivated 
NC is the k-space delocalization (blurring) of electronic states. This effect can be directly 
attributed to the confinement in a finite space by window function wQ^ as discussed in 
the previous chapter. Fig. 12 shows k-space density projection of the lowest unoccupied 
state onto a 1-D line between the r(0,0,0) and r(0,2,2) k-points for different Si:H NC 
diameters. It demonstrates the variation of k-space delocalization and its shape depending 
on the size of the Si NC. We can clearly identify the sinc-like shape envelope corresponding 
to the rectangular step window function (see Fig. 6 for comparison). 

B. SiNCiCHa 

In the next step, we examine an impact of methyl (-CH3) terminating group on the 
atomic and electronic structure of Si NCs. Our fully relaxed SiiCHa NCs expand by ^ 1% 
of lattice constant as a consequence of steric repulsion between the individual -CH3 groups 
(see Fig. 14). The presence of mechanical strain causes inhomogeneities in the atomic Si- 
core structure, which lifts up the degeneracy of molecular states (see Fig. 10b). This effect 
significantly reduces the size of the mini-gaps (see Fig. 10b). Moreover, the electronic states 
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FIG. 14. Histogram of Si-Si bond lengths in Sisas NCs with different passivation. Maps of bond 
length distortions are depicted as insets (a,b,c) in colorscale: (a) Sisas-H, (b) Si538-CH3 and (c) 
Si538-OH. (d) Increase in the mean bond length between -H and -CH3 by 1%. 

are more localized in real space compared to those of a H-passivated NC of a similar size as 
seen from comparison of Fig. 10a, b. Consequently, the tails of projected densities Pi(k') of 
individual electronic states in k-space decay more slowly in k-space than in Si:H NCs of the 
same size (compare again Fig. 10a, b). This means that individual "bands" are not strictly 
localized in their k-momentum. 

We found out that the presence of non-polar methyl groups leads to the localization of 
surface states in energies far from the Si NC band gap similarly as in Si:H NCs (Fig. 9b), 
but the band gap is smaller compared to Si:H NCs (Fig.9d). 

C. SiNCiOH 

According to our DFT simulations. Si NCs capped with hydroxyl (^OH) groups suffer 
from large atomic relaxation in proximity to the surface with respect to Si bulk structure (see 
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Fig. 14c). Nevertheless, atomic structure deep within the Si NC core remains almost intact 
and thus close to the bulk. The large distortion of surface atomic structure is driven by the 
presence of significant charge transfer between Si and electronegative O atoms. The strong 
polarity of the Si-0 bond tends to cause large localization of their MOs near the surface 
in real space (see Fig. 10c) with energies close to band gap edges (Fig. 9c). This picture is 
consistent with previous experimental evidence [30] and theoretical simulations [16-20] of 
oxidized Si NCs. 

What's more, our simulations point out that -OH groups tend to align via hydrogen 
bonds, forming ordered domains at low temperature. Consequently, considerable local elec- 
tric fields are induced across the Si NC as shown in Fig. 13. We computed the electric dipole 
of 120 Debye of the largest -OH passivated NC (3 nm, Sisas), while the electric dipoles of 
the corresponding NCs passivated by -H (0.15 Debye) and -CH3 (0.9 Debye) are negligible. 
The presence of these irregularities in electric field together with the geometrical distortion 
of the surface layer (Fig. 14c) probably cause the irregular real space distribution of HOMOs 
and LUMOs in -OH passivated NCs (Fig. 10c). 

Due to the presence of strongly localized states in the surface layer, the projected band 
structure is significantly blurred (with large tails of projected density y5^(k')), which relaxes 
k-space selection rules of optical transitions and improves the radiative recombination prob- 
ability of indirect transitions. Also, the degeneracy of electronic states is almost removed 
due to the large distortion of atomic structure near the surface. This, however, goes hand in 
hand with higher localization and irregular distribution of MOs in real space, which decrease 
real space overlap of orbitals and thus can limit the radiative recombination probability. 

IV. DISCUSSION 

In this article, we show that band structure can be rigorously defined in nanostructures, 
taking into account several alterations with respect to bulk materials. The first difference is 
that the utilization of the band-structure concept makes sense only for ''larger" nanocrystals, 
with sizes above a certain size limit. For example, for Si NCs this size limit lies between 
1.5 and 2 nm (see Fig. 11), because under this limit the k-space projected MOs do not show 
much of band-like behavior. 

Second, the smaller the nanocrystal, the more pronounced the presence of mini-gaps 
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inside the energy bands, which can hmit the non-radiative relaxation of excitons. For the 
same size nanocrystals, the mini-gaps are wider in more symmetric nanocrystals with weaker 
effect of surface due to the degeneracy of energy levels (see Fig. 10). 

Third, the nanoscale size introduces blurring of the k vector (see Fig. 12 and Fig.4d). Al- 
though this phenomenon is qualitatively easily predictable from the Heisenberg uncertainty 
principle and was already treated e.g. by Hybertsen,[27] our approach allows us to quantify 
the influence of k vector blurring for individual MOs. 

Last, in addition to describing the crystalline core, our computed band structures already 
include the effect of surface states. This is very important, because many types of real-hfe 
semiconductor nanocrystals need to be capped by various surface-terminating groups, which 
profoundly influences their electronic properties. 

V. CONCLUSION 

In this paper, we introduced the general method to map effectively electronic structure of 
aperiodic systems such as NCs from the real space to reciprocal space. This method allows 
us to reassemble the electronic band structure of flnite size systems. We believe that this 
method could provide more insight into the question if the band structure concept can be 
still applied to nanocrystals of different shape, size or chemical composition. In particular, 
we demonstrate that the band structure picture of nm-scale Si NCs can be still adopted 
down to nm size, but with two important consequences of the flnite size of system: (i) 
the discretization (ii) the delocalization of electronic states in the reciprocal space. This 
results also means that efflcient slow red photoluminescence in Si NCs arises from indirect 
X-to-F electron-hole recombination, in agreement with recent experimental evidence [31]. 

We employed this method to investigate the effect of Si NC size and presence of different 
passivation groups including non-polar and polar groups on their band structure. We found, 
that Si NCs caped with methyl group expand by ~1 % due to steric repulsion between 
methyl surface groups. In case of polar hydroxyl group, real space localized states near 
surface form band gap edges, which tends toward rather delocalized k-space states. In other 
words, the band structure preserves, but near the band gap the electronic states becomes 
more blurred. In addition, we found strong alignment of hydroxyl group forming signiflcant 
macroscopic electrostatic dipole moment across Si NC. 
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